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FOREWORD 


This report constitutes a part of the work done on the research 
project entitled "Radiative Transfer Models for Nonhomogeneous Atmosphere 
The work was supported by the NASA Langley Research Center through Grant 
No. NSG-1153. The grant was monitored by Dr. Henry G. Relchle. 
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r;ALUATION OF TRANSMITTANCE OP 
SELECTED INFRARED BANDS 


by 


S. K. Gupta and S. N. Tlwarl 


School of Engineering 
Old Dominion University 
Norfolk, Virginia 23508 


SUMMARY 

Computer programs have been developed for evaluating homogeneous 
path transmittance with line-by-line and quasi-random band model formu- 
lations. Transmit tances for some selecced bands of different gases have 
an obtained using these programs. Selected bands include the CO funda- 
mental band (4.6y), N 2 O 4.5y band, CO 2 4.3y band, U 2 O 6.3y band and the 
CO 2 15y band. Results of theoretical computations are compared with 
available experimental measurements. Significant errors are observed in 
the results obtained from the quasi-random band model formulation, indi- 
cating that it is Inadequate to meet the accuracy requirements for 
atmospheric work. 


1. INTRODUCTION 


The work discussed In the present report developed from our studies 

of the application of radiative transfer models to the detection of gaseous 

pollutants of the atmosphere by remote sensing. It has been shown earlier 

[1]* and Is further affirmed during the present Investigation' that the 

quasi-random band model computation of gas column transmlttances Is far 

more economical than the llne-by-llne model computation. It was considered 

important, therefore, to examine the possibility of using the quasi- 

random band model formulation for transmittance computations as required 

in surface temperature retrieval [2] and other data reduction procedures. 

Transmittance computations were made for several bands of different 

gases using line-by-line and quasi-random band model formulations under 

conditions of temperature and pressure for which experimental measurements 

were available. Theoretical results are compared with the experimental 

results of reference [3]. In most cases, initial computations were made 

with higher resolution and the spectra were then degraded with appropriate 

slit functions for comparison with experimental results. 

Spectral parameters for the lines of the bands considered in the 

present study were obtained from the AFCRL compilation [A, 5] which is 

0 

available at the Langley Research Center. 

The spectral transmittance models are discussed briefly In Sec. 2. 
Detailed comparison of tha transmlttances of the Individual bands is 
presented in Sec. 3. 

*The numbers in brackets indicate references. 


2. SPECTRAL TRANSMITTANCE MODELS 


The spectral transmittance models (llne-by-llne and quasi-random 
band model) and the computation procedures used In the present work have 
been described In detail In [1]. These are discussed here briefly simply 
to direct attention to specific areas of present Interest and to emphasize 
areas of Important differences. 

2.1 Llne-by-Llne Model (Direct Integration) 

Drayson [6] and Kunde and Maguire [7] 'lave described a procedure for 
competing the atmospheric transmittance by using the llne-by-llne model. 

The same procedure, with certain modifications described In [1] , has been 
adapted In the present work. The entire frequency range of Interest Is 
first divided Into a large number of narrow Intervals Ao) . The width Ao) 

Is denoted In >he program by DEL. Each Interval Is then divided Into a 
variable number of sublnter\'als depending upon the number of lines within 
the Interval. Two very narrow subintervals are created on each side of 
a line center. The width of each of these subintervals Is denoted by ALX 
In the program. The transmittance Is computed at four frequency locations 
In each subinterval and Is finally averaged over each interval. 

The monochromatic transmittance at any wavenumber location (*) Is 
given by 

t(o 3) » exp[-<(oj) u] , (2.1) 

where <((i>) Is the total absorption coefficient at the wavenumber u) In 
cm”* - atm"* and u represents the pressure path-length of the absorber 
In cm-atm. 

Total absorption coefficient at any wavenumber (*) consists of 


1 


contributions from all the lines in the vicinity and is computed in two 
parts as 

k(o)) - f (2.2) 

where <^(u) and k^(u) are called the direct and wing contributions 
respectively. Direct contribution originates from lines in very close 
vicinity (on both sides) and is obtained from 

- Z S y /{■’![ (w-w )^ + 1 ) » (2.3) 

n n n n n 

where u) refers to the center of the nth contributing line. In the 
n 

computer program, the range of direct contribution is denoted by DLIM. 

Thus, if 1 (w-'-o^) I £ DLIM , the contribution is called direct and is 
evaluated by using Eq.(2.3). The summation in this equation extends 
over all lines from minus DLIM to plus DLIM. The wing contribution 
arises from lines which are farther from U) than DLIM and is obtained 
from 

- Z Y„/(^(co-w )M . (2.4) 

n n n n 

The wave number range of the wing contribution is denoted by WLIM such that 

DLIM <l(o)-u3 )| < WLIM . Equations (2.3) and (2.4) are based on the 
n — 

assumption that the line shapes are Lorentzian at the tropospheric pressures. 

Different values of the parameters DEL and ALX have been chosen for 
different bands depending on the frequency range of the band and the gas 
pressure. The line-by-line computer program is given the name GASLINE. 

Various symbols used in the program are explained In Appendix A and a 
listing of the program is given in Appendix B-1. 
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2.2 Quasi-Randon Band Modal 

In this formulation, the frequancy range of the band la divided Into 
Intervals (Au>) of equal width which Is denoted In the computer program by 
DEL. The lines within each Interval are distributed Into five Intensity 
decades and average Intensity for each decade Is obtained. Average trans- 
mittance over the Interval Au , due to a single line of Intensity S , 

Is given by [8] 

t(Ao)) - -r- / exp[-S u b(w,o) )ldu , (2.5) 

where £galn u represents the pressure path-length In cm-atm and b(b),u:^) 

Is the shape factor for the line with center at o) . 

0 

£.'.presslon (2.5) can easily be extended, first to Include the contribu- 
tion from all the lines within one Intensity decade and then from the entire 
Interval [8,9]. In this case also, It Is important to take into account 
the contribution from lines In some adjacent intervals. The number of 
adjacent Intervals (on both sides), from which the contribution Is considered 
significant. Is denoted by JD In the program. Different values of DEL 
have been adopted for different bands. The computer program Is given the 
name GASBAND. The symbols used in the program are explained in Appendix A 
and a listing of the program is given In Appendix B-2. 
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3. TRANSMITTANCES OF INDIVIDUAL BANDS 


The bands chosen for the comparison of transmittances In the present 
study are the CO fundamental band and those bands of other minor atmospheric 
constituents which Interfere spectrally with the CO fundamental band. The 
CO 2 ISu band has been chosen because of Its Importance In various atmospheric 
radiative transfer analyses. Transmittance results are presented In Figs. 

3. 1-3. 5 separately for each band. Llne-by-llne results are shown by the 
solid lines, quasi-random band model results by the histograms and the 
experimental results of Burch et al. [3] by the broken lines. Integrated 
absorptance for these bands are presented for comparison in Table 3.1. 

3.1 CO Fundamental (4.6 ) Band 


The frequency range selected for the study of this band (2070-2220 cn“') 
is the same as the range of operation of the gas-filter correlation instrument, 
which is flown from the Langley Research Center. Comparison for this band Is 
shown in Fig. 3.1. A value of D^ . - 1 cm"* was used in the linc-by-llne 
computation. Due to low pressure of the gas in this case (51 mm Hg) , resulting 
in small line widths, a value of ALX 0.003 cm~* was used. For the band 
model computation, a value for DEL ■ 5 cm"* was used. The experimental 
results were obtained with an effective slit-width of 25 cm“* [3]. The 
theoretical results, therefore, were also degraded with a slit function of 
25 cm"* for comparison. It can be seen from Fig. 3.1 that the llne-by-llne 
results are in good agreement with the experimental results while the band 
model results exhibit appreciable differences, particularly in the P and R 
branches of the band. Integrated absorptances for this band, presented in 


Table 3.1, are for slightly greater frequency range than previously considered. 
This was necessary because the experimentally measured absorptance covered a 
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Fig. 3.1 Comparison of transmittance s of CO fundamental band 




greater frequency range. The results presented in Table 3.1 show that 
while the absorptance obtained by the line-by-line model is only five 
per cent lower than the experimental value, it is nit e per cent lower for 
the band model absorptance. 

3.2 N 2 O A.Svi Band 

The frequency range selected for this band is 21A0-2290 cm~‘ and the 
results are shown in Fig. 3.2. For the linc-by-line computation, a value 
of DEL 1 cm~^ was adopted. Since the gas pressure in this case was 
relatively high (644 mm Hg) , ALX ■ 0.02 cm“' was used. For the. band model 
computation, values of DEL - 5 cm”* and JD “ 10 were adopted. The 
experimental results for this band were obtained for an effective slit- 
width .5 cm~* and, therefore, theoretical results were degraded with a 
slit function of 25 cm”* for comparison. It is apparent from Fig. 3.2, that 
the experimental transmittance curve is shifted approximately 6 cra“* to the 
higher frequency side with respect to the theoretical curves. The reason 
for this shift is not clear at this time, h’ith the exception of this, the 
agreement between the line-by-line and experimental results is very good. 

The band model curve again shows appreciably low absorption. It can be seen 
from Table 3.1 that the line-by-line absorptance in less than two per cent 
lower than the experimental value while the difference for the band model 
absorptance is approximately ten per cent. 

3 . 3 COq 4.3u Ban d 

The comparison of experimental and tJieoret lea 1 transmittanecs for this 
band, over the frequency range 2220-2420 cm"*, is shown in Fig. 3.3. Mille 
the linc-by-line results aie obtained for DEL ■ 1 cm”* and /J.X - 0.02 cm"*, 
the band model results are for DEL ■ 5 cm”* and JD • 10. Experimental resul 
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Fig. 3.2 Comparison of transmittances of 4.5 y N 2 O band. 
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Fig. 3.3 Comparison of transmiCtances of 4.3 y C0_ band. 
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shown were obtained for an effective sllt-wldth of approximately 5 cm' 


For purposes of comparison, therefore, the llne-by-llne results were 


slit function while the band model results were 


only averaged over shifted and unshlfted meshes. The agreement between 


the experimental and the llne-by-line results Is seen to be excellent 


The band model curve again exhibits a slightly lower absorption. It can 


be seen from Table 3.1 that the Integrated llne-by-llne absorptance Is 


within 0.5 per cent of the experimental value while the band model 


absorptance is approximately 5.5 per cent lower 


3. A HoO 6.3p Band 


The comparison for this band, which extends over the frequency region 


is shown in Figs. 3.4a - 3.4c. It was found necessary to 


present these results in three figures to maintain clarity. For the fre 


quency range of Fig. 3.4a, the experimental and line-by-llne results are 


shown in Fig. 3.4(a-l) and the quasi-random band results in Fig. 3.4(a-2) 


again to maintain clarity 


The experimental results for this band were obtained with an effective 


slit-width of 6 cm' 


which were 


obtained originally for DEL = 2 cm' 


a 6 cm"^ slit function. To ensure that no significant structure was lest 


attempts were not made to smooth the results any furthe 


and model results 


were obtained for DEL ■ 6 cm' 


shifted and unshifted meshes. No further degradation of these was attempted 


theoretical results in the lower frequency region of the band. In the middle 


and higher frequency regions (Figs. 3.4b and 3.4c), however, the theoretical 
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results exhibit much more structure than the experimental results. Inte- 
grated absorptances for this band, listed In Table 3.1, exhibit excellent 
agreement. 


3.5 CO 2 15y Band 

Comparison of results for this band Is shown In Fig. 3.5. The experi- 
mental results [3] refer to different values of effective sllt-wldth over 
different parts of the band. The llne-by-llne results obtained for DEL > 1 cm" 
and ALX - 0.02 cm~^ were degraded by a 3 cm"^ silt function for the region 
550-650 cm“' while a 5 cm“^ silt function was used for the region 650-800 cm”^. 
The band model results were obtained for DEL <■ 5 cm~^ and JD ■ 10 and were 
averaged over shifted and unshlfted meshes with no further degradation. It 
can be seen that the agreement between the three results Is excellent. 

Table 3.1 shows that the theoretical results agree with the experimental ones 
within one per cent. 

It Is clear from the foregoing comparison that the llne-by-llne 
computation yields results In better agreement with the experimental values 
than the quasi-random band model. In view of the high accuracy required 
for atmospheric work, It Is desirable to use the llne-by-llne model for the 
temperature retrieval work and other data reduction. Also, it Is known that 
the spectra of the various atmospheric constituents under consideration have 
significant correlation among themselves. In the quasi-random band model, 
however, average transmittance over the interval is evaluated as a first step 
and so it is inherently unsuitable for computing transmittances of multi- 
component atmospheres. 




Band Identi- 
fication and 
Freq. Range 

Experimental 
Absorptance 
(Burch et al. 
1962) 

Line-bv-Llne Re«?ults 

Band Mode 

1 Results 

Integrated 
Absorptance 
(cm”' ) 

Percentage 
Diff. with 
Measurement 

Integrated 

Absorptance 

(cm"') 

Percentage 
Diff. with 
measurement 

CO 4.6 M 
(1970-2270 
cm”') 

73.90 

70.32 

-4.84 

67.30 

-8.93 

NoO 4.5 U 
(n40-2290 
cm~' ) 

57.00 

56.01 

-1.74 

51.43 

-9.77 

CO 2 4 . 3 u 
(2220-2420 
cm”') 

73.80 

74.09 

0.39 

69.72 

-5.53 

H 2 O 6.3 u 
(1200-2100 
cm”') 

334.0 

326.4 

-2.28 

323.1 

-3.26 

CO 2 15 u 
(550-800 
cm"' ) 

66.10 

66.57 

0.71 

65.35 

-1.13 






CONCLUSIONS 


Computer programs have been developed for evaluating homogeneous path 


spectral transmittances with line-by-llne and quasi-random band model formu 


latlons. The llne-by-llne program Is given tha nane GASLINE and the band 


model the name GASBAND. By using these programs, transmittances for several 


bands of different gases (for selected pressure and path-length conditions) 


have been computed. The particular bands selected for this study were 


Theoretical 


results were compared with the experimental results of Burch et al. for the 


same set of pressure and path-length conditions. Comparisons of line-by-line 


and experimental results for transmittances indicate good agreement for all 
the bands. The quasi-random band model results for CO (4.6y), N 2 O (A.5y), 
and COj (A.3y) bands, however, exhibit significant differences from the 


experimental results. Comparisons of integrated absorptance for these bands 


indicate Che general trend observed for the transmittances 
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APPENDIX A 

EXPLANATION O'.’ SYMBOLS USED IN THE 
COMPUTER PROGRAMS 


A-1. Synbolt Used in Program GASLINE 

AC Monochromatic absorption coefficient at the 

specified frequency. 

ACB, ACE I Total wing contribution to the absorption coefficient 

& ACM I at the interval boundaries and center. 

AL Half-widths of the Individual lines. 

ALX Width of a narrow sub-interval near the line center. 

DEL Width of an Interval. 

DLIM Width of the region from which the direct contribution 

Is obtained. 

EL Energy of the lower states for the Individual lines. 

FR Frequencies of the individual lines. 

FRB Frequencies at the Interval boundaries. 

FRC Frequencies at the Interval centers. 

FRE Frequencies of the lines falling within one Interval. 

FRC Frequencies of the individual mesh points within one 

Interval. 

FRL Lower frequency boundary of the band. 

FRS Frequencies at the sub-interval boundaries. 

Fr'? Upper frequency boundary of the band. 

KR Number of intervals in the band. 

L£ Total number of lines in the spectrum. 
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WLl, WL2 


WLIM 


XLl, XL2 





Number of lines Ic one intervel. 

Number of sub-interval boundaries. 

Number of sub-intervals in one Interval. 

Pressure path-length of rhe absorber. 

Pressure referring to NIP (760 mi> Hg). 

Effective pressure of the absorber. 

Exponent to account for the tenperaturc-dependenv'.e 
of the rotational partition function. 

Intensities of the individual lines. 

Slopes for the wing contribution between 
interval boundaries and the center. 

Temperature of the absorber (*K). 

Reference temperature for the molecular 
spectral parameters (*K). 

Temperature referring to NTP (273*K). 

Monochromatic transmittance at the specified frequency. 

Transmittance averaged over each interval. 

Transmittance averaged over wider intervals to account 
for the slit functions. 

Factor accounting for the temperature-dependence of 
the vibrational partition function. 

Weighting factors for the Gauss-LegenJre quadrature 
formula. 

Width of the region from which the wing contribution 
is obtained. 

Abscissa values for the Gauss-Lcgendre quadrature formula. 


RfTROmTIBII IT.' (' ' : 

UUICINA.. ■ I 





! 


f 


I 


4 

i 




i 


C 


t 


! 







23 





A-2. Symbols Used In Program GASBAND* 


Average half-width of the lines of the molecule. 


Average of the intensities of the lines within one 
intensity decade. 


Upper intensity limit for individual intensity decades. 


Intensity of the strongest line within one Interval. 


Index which lets the transmittances be computed for 


the unshifted and shifted meshes. 


Number of adjacent Intervals on both sides of an 


interval from which the contribution is taken into 


account. 


Serial number of the first line in an interval. 


Serial number of the last line in an Interval. 


Number of lines within one decade of an Interval 


Value of the integral for one line - direct or wing 


contribution. 


Sum of the intensities of all the lines within one 


decade of an interval. 


Weighting factors used in the nunorlcal integration of 


direct contribution. 


Weighting factors used in the numerical integration of 
wing contribution. 


Array of transmittance values containing results from 
unshifted- and shifted-mesh computations. 


Contribution to the total transmittance from lines 


within one interval. 


*As far as possible, same symbols have been used in the Program 
GASBAND as in GASLINE. Only the symbols which are defined diff- 


erently for GASBiiND are explained here. 
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APPENDIX B 


LINE-BY-LINE AND QUASI-RANDOM BAND 
MODEL COMPUTER PROGRAMS 


B-1. Program GASLINE 


read 

read 

read 

READ 

read 


10 
1 1 
12 
13 



program gasline ( input •output ) 
integer x*y 

DIMENSION FR (4590 ) »S I (4590 ) »EL( A590 1 »AL (4590 ) •P’RB (451 ) .FRC(450> • 
/FRE( 100 ) .FRS( 126) .FRG(500) . AC ( 500 ) % TR ( 500 ) * TRA ( 450 ) ♦TRAV(450) 


FRL . FRU • DEL • OL I M , WL I M 
PREC , PNTP • TEMC • TEMR « TNTP 
XLl .XL 2 , WLl •WL 2 
PL»VP. ALX»RPtLE 

(FR(X) tSl (X)«EL(X) .AL(X)*X = 1 tLE) 


10 
10« 

1 1 « 

12 . 

13. 

FORMAT ( 10F8.2 ) 
format (8F10.6) 

FORMAT (4F8. 4 . 18 ) 

FORMAT (FI 1.3.E13.4.F11.3.F5.3) 

PI =3. 14 159 

LINE INTENSITIES AND HALF-WIDTHS ARE CONVERTED 
AMBIENT CONDITIONS OF TEMPERATURE AND PRESSURE 
CST= ( SORT ( TEMR/TEMC ) ) *PREC/PNTP 
PART=VP*{ (TEMR/TEMC )**RP) *2. 69E+ 19 
FACT= 1 .439* ( TEMC-TEMR)/ (TEMC*TEMR) 

DO lOl X=1.LE 
AL(X )=CST*Al(X) 

101 SI (X)=SI (X)*PART*EXP(EL(X)*FACT) 

ENTIRE frequency RANGE IS DIVIDED INTO INTERVALS 
AT interval BOUNDARIES AND CENTERS ARE DEFINED 


from 


FREQUENCIES 


DELA=0.5*DEL 

RK= (FRU-FRL )/DEL+0. I 

KR = RK 

FR3( 1 )=FRL 
DO 102 K=1 ,KR 
FR3(K+1 )=FRB(K)+DEL 
102 FRC (K ) =FR3 (K )+DELA 
ALYr:2. *ALX 

C LOOP BELOW GOES TO THE END OF THE PROGRAM - COVERS ALL STEPS 
Y = 0 

DO 103 K= 1 ,KR 

C DETERMINES THE NUMBER OF LINES (MR) IN THE PARTICULAR INTERVAL 
M = 0 
MP = M 

106 IF (Y-LE) 104, I 00 .100 
104 Y=Y+1 

I F (FR{ V )-FR3 (K ) ) 106.107.107 

107 IF(FR(Y )-FRB(K+l ) ) 108.105. 105 

ICO M=M+1 

FRE (M ) sFR ( Y ) 
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MP=M 
GO TO 106 
105 Y=Y-1 
100 CONTINUE 

DIVIDES THE INTERVAL INTO SUB- 1 NTERVALS AND UP-DATES THE NUMBER OF 
SUB- INTERVAL BOUNDARIES NP 
N^l 

FRS(N)--FRB(K) 

IF (MP) 109.109«110 
no CONTINUE 

DO 111 M=1*MP 
DIF=FRE(M)-FRS(N) 

IF (OIF-ALX) 112.112*113 

113 IF (DIF-ALY) 114.114.115 
I 12 FRS(N+1 )=FRE(M) 

N = N+1 
GO TO 1 16 

114 FRS(N+1 )=FRE(M)-ALX 
FRS(N+2)=FRE(M) 

N = N+2 
GO TO 116 

115 FRS(N+1 )=FRE(M)-ALY 
FRS { N-»-2 ) =FRE { M ) -ALX 
FR5(N+3)=FRE(M) 

N = N+3 

116 IF (M-MP) 117. I 18.118 

117 D1F=FRE(M+1 )-FRS{N) 

IF (DIF-ALX) 119.119.120 

120 IF (DIF-ALY) 121. 121.122 
1 19 GO TO 1 1 1 

121 FRS(N+1 )=FRE(M)ValX 
N = N+1 

GO TO 111 

122 FRS(N+I )=FRE (M )+ALX 
FRS ( N+2 ) =FRE ( M ) •♦■ALY 
N = N+2 

111 CONTINUE 

118 DIF=FRB(K+1 )-FRE(M) 

IF (DIF-ALX) 123,123.124 

124 IF (DIF-ALY) 125.125,126 

123 FR5(NH )=FR3(K+1 ) 

NP = N-H 

GO TO 127 

125 FRS(N-H )=FRE(M)+ALX 


{ 
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FRS(N+2 )=FR0(<-M ) 

NP=N+2 
GO TO 127 

126 FRSlN+1 )=FRE IM )+ALX 
FRS ( N+2 ) =FRE { M ) +ALY 
FRS(N+3 )=FR3(K+1 ) 

NP=N+3 

GO TO 127 

109 FRS(N+1 )=FRB(K-M ) 

NP=N+1 

127 NQ=NP-1 

C generates the mesh by COMPUTING FOUR FREQUENCY LOCATIONS IN EACH 

c SUB- interval 

DO 128 N=1 «NQ 

VAR=0.5* (FRS(N+l )-FRS(N) ) 

C0N=0.5*(FRS(N+I )+FRS(N>) 

I=4*(N-1 ) + l 
FRG( I ) =CON-VAR*XLI 
FRG( I + l )=C0N-VAR*-XL2 
FRG ( I +2 ) =C0N+VAR*XL2 

128 FRG( 1+3 )=CON+VAR»XLI 
IG=4*NQ 

C INITIALIZES ABSORPTION COEFFICIENTS AT THE MESH POINTS AND THE 
C boundaries and centers of THE INTERVAL 
ACB=0. 

ACM=0* 

ACE=0. 

DO 129 1=1. IG 

129 AC ( I )=0* 

C STARTS EVALUATION OF THE ABSORPTION COEFFICIENT 
DO 130 X=1.LE 
DIF=ABS(FR(X)-FRC(K) ) 

IF (DIF-WLIM) 13 1.1 31 .130 

131 IF (OIF-DLIM) 132.132.133 

C evaluates DIRECT CONTRIBUTION AT THE MESH POINTS 

132 DO 134 1=1 . IG 

FD=FR ( X )-FRG( I ) 

DEN = Pl>(FD*FD+AL(X)»ALfX) ) 

134 AC( I )=AC( I )+SI (X)*ALIX)/DEN 
GO TO 130 

C evaluates wing CONTRIBUTION AT INTERVAL BOUNDARIES AND CENTER 

133 FB=FR (X ) -FRB (K ) 

FM = FR(X )-FRC(K) 

FE=r-R(X )-FR3(K+l ) 







PNUM=S1 (X)*AL(X) 

ACB= AC0+PNUM/ ( P I *F0*FB ) 

ACM= ACM+PNUM/ ( P I *FM#FM ) 

ACE=ACE+PNUM/ (PI *FE*FE ) 

130 CONTINUE 

C COMPUTES SLOPES FOR W I NG CONTRIBUTION BETWEEN INTERVAL 
C BOUNDARIES AND CENTER 

SLl = (ACM-AC3)/(rRC(K)-FRB(K) ) 

SL2= (ACE- ACM)/ (FR3(K+1 )-FRC(K) ) 

C EVALUATES TOTAL ABSORPTION COEFFICIENT AT ALL MESH POINTS 
DO 136 1=1 tIG 
DIF=FRG( I )-FR3(K) 

IF (DIF-DELA) I37tl3a»138 

137 AC( I )=AC( I )+ACO+SLl*OlF 
GO TO 136 

138 AC( I )=AC( 1 )+ACM+SL2* (DIF-DELA ) 

136 CONTINUE 

C COMPUTES MONOCHROMATIC TRANSM I TT ANCES AT ALL MESH POINTS 



DO 1 4 I I = I . I G j 

00=PL*AC(I) i 

IF (0D.GT.675.) GO TO 153 

TR ( I )=EXP(-OD) ^ 

GO TO 141 
153 TR ; I )=0. 

141 CONTINUE 

C EVALUATES TRANSM I TTANCES AVERAGED OVER INTERVAL 
TRA(K)=0. 

DO 144 N=1,NQ 
VAR=0.5*(FRS(N+l )-FRS(N) ) 

I =4* (N-I )^-l 
SUMl =TR ( I )+TR( 1+3 ) 

SUV2=TR( I+l )+TR( 1+2) 

SUM=WL1 +SUMI +WL2*SUM2 
144 TRA (K ) =TRA (< )+SUM*VAR/DEL 
103 CONTINUE 

C SIMULATES SLIT FUNCTION TREATMENT BY TAKING 5-POlNT SLIDING AVERAGES 
i, KN=KR-2 

DO 151 K = 3.KN 

f 151 TRAV (K ) = (TRA (K-2 ) +TRA (K-1 ) +TRA (K )+TRA (K+1 )+TRA (<+2 ) )/5. 

PRINT 60» ( TRA ( K ) »K= 1 »KR ) 

PRINT 60« (TRAV(K ) .K=3.KN) 

, 60 FORMAT! 1H1/////(10F10«5»/) ) 

STOP 


END 
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REPRODUCIBn.ITY OF THE 
ORIGINAL PAGE IS PO(V! 


2. Prot>ram GASBAND 


program GAS0AND ( input I output ) 
integer XiW 

dimension FR(4590 ) «SI (4590 ) .EL (4590 ) .FRB( 152 ) »FRC( I51).LIB(151). y 
/LIE(151).MP (151). BIGI(151).3IG(6. 151). SS 1(5.151). NS 1(5. 151)* 

/A VS I (5.151) .TRG( 151 ) .TRA(301 ) » X 1 ( 26 ) . T 1 ( 26 ) . X2 ( 2 1 ) . T2 ( 2 1 ) 

READ 10. FRL.FRU.DEL .°NTP.PREC.TNTP.TEMC.TEMR j 

read 11. LE.JO 

read 12. PL. VP. RP 

read 12. (XI ('*') iN = l .26) ; 

READ 12. (T1 (W) .A’ = l .26) ’ 

READ 12 . (X2 (W ) .W=l .21 ) I 

read 12 . (T2 (W ) .W=1 .21 ) 

READ 13. (FR (X ) .S I (X ) .EL (X ) .X=l .LE ) 

I 0 FORMAT ( lOFB. 1 ) 

II FORMAT (1615) 

12 FORMAT { 1 OFB. 4 ) 

1 3 FORMAT (2 (FI 1.3.EI3.4,F11.3.5X>) 

CONVERTS THE LINE INTENSITIES AND THE AVERAGED LINE-WIDTHS FROM 
REFERENCE TO AMBIENT CONDITIONS OF TEMPERATURE AND PRESSURE 
CST= (SORT (TEMR/TEMC ) )>?REC/PNTP 
PART = VP* ( ( TEMR/TEV.C ) **RP)*2,69E+\9 
FACT= 1 .439^ (TEMC-TEMR)/(TEMC»T£MR) 

DO 101 X=1.LE 

101 SI(X)=SI(X) -»PART*EXP (EL(X )*FACT ) 

ALA=AL»CST 

INITIALIZES THE FREQUENCY INTERVAL BOUNDARIES FOR THE UNSHIFTED MESH 
DELA=0.5*DEL 
RK= (FRU-FRL )/DEL+0. 1 
KR = RK 

FRB( 1 )=FRL 
IR=0 

STARTS the COMPUTATION FOR THE UNSHIFTED MESH 
GO TO 121 
122 KR=KR-1 

FRO( 1 )=FRL-DELA 

DEFINES THE FREQUENCIES AT THE INTERVAL BOUNDARIES AND CENTERS 
121 DO 1 00 K=1 .<R 

FR3(K+1 )=FRB <K )+DEL 
100 FRC(<)=FRa(K )+D£LA 

DETERMINES THE NUMBER OF LINES IN EACH FREQUENCY INTERVAL 
X = 0 

DO 1 02 K-l .KR 
M=0 

104 IF (X.GE.LF) GO TO 103 
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X=X4-1 

IF (FR(X ) .LT.FQB(<) ) GO TO 104 
IF (FR/X) .GE.FR3(K+1 ) ) GO TO 103 
M = M+1 
GO TO 104 

105 IF (K.GE.KR) GO TO 103 
X = x-1 

103 MP(K)=M 

102 CONTINUE 

ASSIGNS NUMBERS TO THE FIRST AND LAST LINE OF EACH INTERVAL 
X = 0 

DO 106 K=1 •KR 

IF (MP(K).EQ*0) GO TO 106 

LIB(K)=X+1 

LIE(K)=LIB(K )+MP(K)-l 
X=L1E(K ) 

106 CONTINUE 

FINDS OUT THE HIGHEST INTENSITY IN THE INTERVAL* GENERATES THE 
intensity decades* DISTRIBUTES LINES INTO INTENSITY DECADES* 
DETERMINES THE NUMBER OF LINES IN EACH DECADE AND THEN COMPUTES 
average intensity p-OR each DECADE 
DO 1 07 K=1 ,<R 
JB=LIB (K) 

JE=LIE(K) 

BIGI (K)=SI ( JB) 

JC=J0+1 

DO 1 08 J = JC*JE 

IF (BIGI (K) .GE.SI (J) > GO TO 108 
BIGI (K)=Sr ( J) 

103 CONTINUE 

DO 109 1=1*6 
IX=-I+1 

109 BIG( I *.<)=BIGI (K)*l 0.**IX 
00 110 1=1*5 
N = C 

SSI(1*K)=0. 

DO 111 J=JB*JE 

IF (3 I ( U ) *GT.r3 I G ( I .< ) ) GO TO 111 
1F(SI ( J) .LE.BIGI 1+1 *K) ) GO TO 111 
N = N+1 

SSI ( I *<)=SSI ( I *0+51 (J) 

NS I ( I *0 =N 
11! CONTINUE 

IF (NS I ( I *K ) *GE . 1 ) 


GO TO 110 
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PFTKOPOT'r 11’'^ • 

OBICINA' I 

NSI ( I «K )=1 

no AVSl ( I *K )=SSt ( I tO/NSl ( I «K ) 

107 CONTINUE 

PI=3.1A159 
RH0=ALA/0ELA 

STARTS the computation OF TRANSM I TTANCES FOR EACH I 
DO 112 K=l »KR 
TRGIK )=1 • 

DO 113 J=lfKR 
TRD=1 • 

JA= I A0S( J-K) 

IF (JA.GT.JD) GO TO 113 
ZI =FRC (K )-FRC( J) 
epsi=zi/dela 
DO 115 I =l »5 
NSJ=NS I ( I » J ) 

PNUM=RHO*RHO*AVS I ( I t J ) »PL/ ( P I *ALA ) 

RES=0* 

IF ( J.NE.K ) GO TO 116 

EVALUATES THE INTEGRAL FOR DIRECT CONTRIBUTION 
DO 117 W=1 .26 

YY=PNUM/(Xl (W)*X1 (W)+RHO»RHO) 

IF (YY.GT.675.) GO TO 119 
Y=EXP(-YY ) 

GO TO 117 

1 19 Y = 0 ^ 

1 1 7 RES=RES+Y»T1 (W ) 

GO TO 115 

evaluates the integral for WING CONTRIBUTION 
116 DO 1 la W=1 .? 1 

YY = PNUM/( (E?SI-X2(W) ) * ( EPS I -X2 ( W ) ) ) 

IF (YY.GT.675.) GO TO 120 
Y-EXP(-YY) 

GO TO 118 

120 Y=0. 

118 RES = RES + Y*T2 { W ) 

RES=RES/6. 

115 TRD=TRD»RES**NSJ 
113 TRG(K )-TRG(K )*TRD 
1 12 CONTINUE 
IR-IR-M 

IF { IR.GT . 1 ) GO TO 123 
DO 1 1 A K= I .KR 
llA TRA(2iK-l ) = TRG(K) 

GOES ro 122 A^JO STARTS COMPUTATION FOR THE SHIFTED MESH 
GO TO 122 
123 DO l?A K=1 .KR 
12A TRA ( ?*< ) = TRG (<) 

AVERAGES RESULTS OVER SHir-ED AND UNSHIFTED MESHES 
DO 125 K- I ,<R 

125 TrRG(K ) = (TRA { 2»K-1 ) f TRA ( 2»K ) + TRA (2»K+1 ) )/3. 

PRINT 60. ( TRG (<).<- 1 . KR ) 

60 FORMAT ( IHI//// ( 10ri2.5.//) ) 

STOP 
END 
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